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Preface 


The  work  described  herein  was  conducted  at  the  U.S.  Army  Engineer 
Waterways  Experiment  Station  (WES),  Coastal  and  Hydraulics  Laboratory 
(CHL)  as  part  of  the  Coastal  Inlets  Research  Program  (CIRP)  under  the  work 
unit  Modeling  Waves  at  Inlets.  Overall  program  management  for  CIRP  is 
directed  by  the  Hydraulic  Design  Section  of  the  Headquarters,  U.S.  Army  Corps 
of  Engineers  (HQUSACE).  Program  Monitors  for  the  CIRP  are  Messrs.  Barry 
W.  Holliday,  John  Bianco,  and  Charles  B.  Chesnutt,  HQUSACE.  The  Program 
Manager  is  Mr.  E.  Clark  McNair,  CHL,  and  CIRP  Technical  Manager  is 
Dr.  Nicholas  C.  Kraus,  CHL. 

The  purpose  of  the  Modeling  Waves  at  Inlets  work  unit  is  to  provide  field 
tools  for  quantifying  wave  processes  at  coastal  inlets.  Wave  information  is 
required  for  almost  all  engineering  studies  near  inlets  to  estimate  channel 
shoaling  and  migration,  shoreline  change,  scour,  forces  on  structures,  and 
navigation  safety.  This  report  documents  a  “workhorse”  numerical  model  for 
estimating  nearshore  wave  growth  and  transformation.  The  purpose  of  the  report 
is  to  transfer  this  technology  to  field  users,  through  guidance  on  model 
application. 

This  study  was  conducted  by  CHL  personnel  under  the  general  direction  of 
Dr.  James  R.  Houston,  Director,  and  Mr.  Charles  C.  Calhoun,  Jr.  (retired). 
Assistant  Director.  Direct  guidance  was  provided  by  Messrs.  Thomas  W. 
Richardson,  Chief,  Coastal  Sediments  and  Engineering  Division,  and 
Mr.  Bruce  A.  Ebersole,  Chief,  Coastal  Processes  Branch  (CPB).  This  report  was 
prepared  by  Drs.  Jane  McKee  Smith,  CPB,  and  Donald  T.  Resio,  Senior 
Research  Scientist,  CHL,  and  Dr.  Alan  K.  Zundel,  Engineering  Computer 
Graphics  Laboratoiy,  Brigham  Young  University.  Assistance  and  review  were 
provided  by  Mr.  S.  Jarrell  Smith,  CPB,  Dr.  Joon  P.  Rhee,  Prototype 
Measurement  and  Analysis  Branch,  and  Ms.  Lori  Hadley,  Coastal 
Hydrodynamics  Branch. 
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1  Introduction 


Estimating  nearshore  wind-wave  growth  and  transformation  is  a  critical 
component  of  most  coastal  engineering  projects,  e.g.,  predicting  bathymetric  and 
shoreline  change,  estimating  navigation  channel  shoaling  and  migration, 
designing  or  repairing  coastal  structures,  assessing  navigation  conditions,  and 
evaluating  natural  evolution  of  coastal  inlets  or  beaches  versus  consequences  of 
engineering  actions.  Nearshore  wave  propagation  is  influenced  by  complex 
bathymetry  (including  shoals  and  navigation  channels);  tide-,  wind-,  and  wave¬ 
generated  currents;  tide-  and  surge-induced  water-level  variation;  and  coastal 
structures.  Use  of  numerical  wave  models  has  become  widespread  to  represent 
wave  transformation  primarily  because  of  their  increasing  sophistication  and 
economy  of  application  relative  to  the  large  expense  of  field  measurements  or 
physical  model  studies. 

This  report  describes  the  application  of  the  steady-state  spectral  wave  model, 
STWAVE.  The  purpose  of  ST  WAVE  is  to  provide  an  easy-to-apply,  flexible, 
and  robust  model  for  nearshore  wind-wave  growth  and  propagation.  Recent 
upgrades  to  the  model  include  wave-current  interaction  and  steepness-induced 
wave  breaking.  STWAVE  has  also  been  incorporated  into  the  Surface-Water 
Modeling  System  (SMS),  which  provides  a  user  interface  and  supporting 
software  for  grid  generation,  interpolation  of  current  fields,  generation  of  input 
spectra,  and  visualization  of  model  output. 

This  report  describes  procedures  for  using  STWAVE.  In  Chapter  2,  an 
overview  of  the  model  governing  equation  and  the  numerical  discretization  is 
presented.  In  Chapter  3,  the  necessary  input  files  are  described.  Guidelines  for 
selecting  parameter  values  are  also  given.  In  Chapter  4,  the  content  and  format 
of  the  STWAVE  output  files  are  described.  Chapter  5  gives  a  tutorial  on  SMS 
application.  Chapter  6  provides  two  example  applications  of  the  model  for  wave 
growth  and  transformation  over  simple  bathymetry  and  wave  transformation  over 
complex  bathymetry  with  wave-current  interaction.  Sample  input  files  are  given 
in  Appendixes  A-D. 
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2  Governing  Equations  and 
Numerical  Discretization 


This  chapter  gives  an  overview  of  the  phase-averaged  spectral  wave  model 
STWAVE  (STeadv-state  spectral  WAVE  model)  (Resio  1987, 1988a,  1988b; 
Davis  1992).  STWAVE  is  a  steady-state  finite-difference  model  based  on  the 
wave  action  balance  equation.  This  report  describes  STWAVE  version  2.0, 
which  includes  the  addition  of  wave-current  interaction. 


Model  Capabilities 

The  purpose  of  applying  nearshore  wave  transformation  models  is  to  describe 
quantitatively  die  change  in  wave  parameters  (wave  height,  period,  direction,  and 
spectral  sh^)  between  the  offshore  and  the  nearshore  (typically  depths  of  20  m 
or  less).  In  relatively  deep  water,  the  wave  field  is  fairly  homogeneous  on  the 
scale  of  kilometers;  but  in  the  nearshore,  where  waves  are  strongly  influenced  by 
variations  in  bathymetry,  water  level,  and  current,  wave  parameters  may  vary  sig¬ 
nificantly  on  the  scale  of  tens  of  meters.  Offshore  wave  information  is  typically 
available  from  a  wave  gauge  or  a  global-  or  regional-scale  wave  hindcast  or  fore¬ 
cast  Nearshore  wave  information  is  required  for  the  design  of  almost  all  coastal 
engineering  projects.  Waves  drive  sediment  transport  and  nearshore  currents, 
induce  wave  setup  and  runup,  excite  harbor  oscillations,  or  impact  coastal  struc¬ 
tures.  The  longshore  and  cross-shore  gradients  in  wave  height  and  direction  can 
be  as  important  as  the  magnitude  of  these  parameters  for  some  coastal  design 
problems. 

STWAVE  simulates  depth-induced  wave  refraction  and  shoaling,  current- 
induced  refinction  and  shoaling,  depth-  and  steepness-induced  wave  breaking, 
diffraction,  wave  growth  because  of  wind  input,  and  wave-wave  interaction  and 
white  capping  that  redistribute  and  dissipate  energy  in  a  growing  wave  field. 

A  wave  spectrum  is  a  statistical  representation  of  a  wave  field.  Concepmally, 
a  spectrum  is  a  linear  superposition  of  monochromatic  waves.  A  spectrum 
describes  the  distribution  of  wave  energy  as  a  function  of  frequency  (one¬ 
dimensional  spectrum)  or  frequency  and  direction  (two-dimensional  spectrum). 


Chapter  2  Governing  Equations  and  Numerical  Discretization 


An  example  of  a  one-dimensional  wave  spectram  is  given  in  Figure  1.  The  peak 
period  of  the  spectrum  is  the  reciprocal  of  the  frequency  of  the  peak  of  the  spec¬ 
trum.  The  wave  height  (significant  or  zero-moment  wave  height)  is  equal  to  four 
times  the  area  under  the  spectrum.  For  the  example  spectrum  given  in  Figure  1, 
the  peak  frequency  is  0.105  Hz,  the  peak  period  is  9.5  sec,  and  the  wave  height  is 
2.8  m.  STWAVE  is  based  on  the  assumption  that  the  relative  phases  of  the  spec¬ 
tral  components  are  random,  and  thus  phase  information  is  not  tracked  (i.e.,  it  is  a 
phase-averaged  model).  In  practical  applications,  wave  phase  information 
throughout  a  model  domain  is  rarely  known  accurately  enough  to  initiate  a  phase¬ 
resolving  model.  Typically,  wave  phase  information  is  only  required  to  resolve 
wave-height  variations  near  coastal  structures  for  detailed,  near-field  reflection 
and  diffraction  patterns.  Thus,  for  these  situations,  a  phase-resolving  model 
should  be  applied. 


Rgure  1 .  Sample  one-dimensional  wave  spectrum 
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Model  Assumptions 

The  assumptions  made  in  STWAVE  version  2.0  are  as  follows; 

a.  Mild  bottom  slope  and  negligible  wave  reflection.  STWAVE  is  a  half¬ 
plane  model,  meaning  that  wave  energy  can  propagate  only  from  the  off¬ 
shore  toward  the  nearshore  (±87.5  deg  from  the  x  axis  of  the  grid,  which  is 
typically  the  approximate  shore-normal  direction).  Waves  reflected  from 
the  shoreline  or  from  steep  bottom  features  travel  in  directions  outside  this 
half  plane  and  thus  are  neglected.  Forward-scattered  waves,  e.g.,  waves 
reflected  off  a  structure  but  traveling  in  the  +x  direction,  are  also 
neglected. 

b.  Spatially  homogeneous  offshore  wave  conditions.  The  variation  in  the 
wave  spectrum  along  the  offshore  boundary  of  a  modeling  domain  is 
rarely  known,  and  for  domains  on  the  order  of  tens  of  kilometers,  is 
expected  to  be  small.  Thus,  the  input  spectrum  in  STWAVE  is  constant 
along  the  offshore  boundary.  Future  versions  of  the  model  may  allow 
variable  input 

c.  Steady-state  waves,  currents,  and  winds.  STWAVE  is  formulated  as  a 
steady-state  model.  A  steady-state  formulation  reduces  computation  time 
and  is  appropriate  for  wave  conditions  that  vary  more  slowly  than  the  time 
it  takes  for  waves  to  transit  the  conq)utational  grid.  For  wave  generation, 
the  steady-state  assumption  means  that  the  winds  have  remained  steady 
sufiBciently  long  for  the  waves  to  attain  fetch-limited  or  fully  developed 
conditions  (waves  are  not  limited  by  the  duration  of  the  winds). 

d.  Unear  refraction  and  shoaling.  STWAVE  incorporates  only  linear  wave 
refraction  and  shoaling,  thus  does  not  represent  wave  asymmetry.  Model 
accuracy  is  therefore  reduced  at  large  Ursell  numbers  (wave  heights  are 
underestimated). 

e.  Depth-uniform  current.  The  wave-current  interaction  in  the  model  is 
based  on  a  current  that  is  constant  through  the  water  column.  If  strong 
vertical  gradients  in  current  occur,  their  modification  of  refraction  and 
shoaling  is  not  represented  in  the  model.  For  most  applications,  three- 
dimensional  current  fields  are  not  available. 

/.  Bottom  friction  is  neglected.  The  significance  of  bottom  friction  on  wave 
dissipation  has  been  a  topic  of  debate  in  wave  modeling  literature.  Bot¬ 
tom  fiiction  has  often  been  applied  as  a  tuning  coefficient  to  bring  model 
results  into  ahgmnent  with  measurements.  Although  bottom  fiiction  is 
easy  to  apply  in  a  wave  model,  determining  the  proper  fiiction  coefficients 
is  chfficult.  Also,  propagation  distances  in  a  nearshore  model  are  rela¬ 
tively  short  (tens  of  kilometers),  so  that  the  cumulative  bottom  fiiction 
dissipation  is  small.  For  these  reasons,  bottom  fiiction  is  neglected  in 
STWAVE. 
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Ongoing  research  will  enhance  present  model  capabilities  and  eliminate  some 
model  assumptions.  The  following  sections  describe  wave  propagation  and 
source/sink  terms  in  STWAVE  version  2.0. 


Governing  Equations 

Interaction  of  waves  with  currents  is  considered  in  a  reference  frame  moving 
with  the  current  Wave  parameters  in  this  frame  are  denoted  with  the  subscript  r, 
for  being  “relative”  to  the  current,  and  parameters  in  the  nonmoving  reference 
frame  are  subscripted  a,  for  “absolute.”  The  wave  dispersion  relationship  is 
given  in  the  moving  reference  frame  as  (Jonsson  1990  and  others) 

ti)^  =  gk  tanh  kd  (1) 

where 

0)  =  angular  frequency 
g  =  gravitational  acceleration 
k  =  wave  number 
d  =  water  depth 

In  the  absolute  frame  of  reference,  the  dispersion  equation  is 

(0  =0)  +  kU  cos(6  -  a)  (2) 

where 

U  =  current  speed 

6=  direction  of  current  relative  to  a  reference  frame  (the  x  axis,  here) 

a  =  wave  orthogonal  direction  (normal  to  the  wave  crest)  (Figure  2) 

The  wave  number  is  solved  by  substituting  Equation  1  into  Equation  2  and  itera¬ 
tively  solving  for  k.  The  wave  number  and  wavelength  (L=(25)/fe)  are  the  same 
m  both  reference  frames. 

Solutions  for  refraction  and  shoaling  also  require  wave  celerities,  C,  and  group 
celerities,  Cg,  in  both  reference  frames.  In  the  reference  frame  relative  to  the 
current. 
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c,.  =  0.5C, 


2kd 

sinh  2kd 


(3) 

(4) 


The  direction  of  both  the  relative  celerity  and  group  celerity  is  a,  the  wave  ortho¬ 
gonal  direction.  In  the  absolute  reference  frame, 

+  U  cos  (6  -  a)  (5) 

(C„),  =  (cp,  ^  m,  <6) 


where  i  subscript  is  tensor  notation  for  the  x  and  y  components.  The  direction  of 
the  absolute  celerity  is  also  in  the  wave  orthogonal  direction.  The  absolute 
group  celerity  defines  the  direction  of  the  wave  ray,  so  the  wave  ray  direction 
(Figure  2)  is  defined  as 
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The  distinction  between  the  wave  orthogonal  (direction  perpendicular  to  the 
wave  crest)  and  the  wave  ray  (direction  of  energy  propagation)  is  important  in  , 
describ-ing  wave-current  interaction.  Without  currents,  the  wave  rays  and 
orthogonals  are  the  same,  but  with  currents,  the  wave  energy  moves  along  the 
rays  whereas  the  wave  direction  is  defined  by  the  orthogonals. 

The  wave  orthogonal  direction  for  steady-state  conditions  is  given  by  (Mei 
1989;  Jonsson  1990) 

C  Da  ^  _  CJc  Dd 

DR  sinh  2kd  Dn  k  Dn 


where  £>  is  a  derivative,  is  a  coordinate  in  the  direction  of  the  wave  ray,  and  n 
is  a  coordinate  normal  to  the  wave  orthogonal. 

The  governing  equation  for  steady-state  conservation  of  spectral  wave  action 
along  a  wave  ray  is  given  by  (Jonsson  1990): 


(9) 


where 

E  =  wave  energy  density  spectrum  (which  is  a  function  absolute  angular 
frequency  and  direction  0) 

S  =  energy  source  and  sink  terms 


Refraction  and  shoaling 

Refraction  and  shoaling  are  implemented  in  STWAVE  by  applying  the 
conser-vation  of  wave  action  along  backward  traced  wave  rays.  Rays  are  traced 
in  a  piecewise  manner,  from  one  grid  column  to  the  next.  The  two-dimensional 
wave  spectra  are  set  as  input  along  the  first  grid  column  (the  offshore  boundary). 
For  a  point  on  the  second  grid  column,  the  spectrum  is  calculated  by  back  tracing 
a  ray  for  each  frequency  and  direction  component  of  the  spectrum.  The  ray 
direction,  a  is  determined  by  Equation  7.  Only  ray  directions  propagating 
toward  the  shore  (-87.5  to  +87.5  deg)  are  included.  Energy  propagating  toward 
the  offshore  is  neglected. 
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The  wave  ray  is  traced  back  to  the  previous  grid  column,  and  the  length  of  the 
ray  segment  DR  is  calculated.  Derivatives  of  depth  and  current  components  nor¬ 
mal  to  the  wave  orthogonal  are  estimated  (based  on  the  orthogonal  direction  at 
column  2)  and  substituted  into  Equation  8  to  calculate  the  wave  orthogonal 
direc-tion  at  column  1 .  Then,  the  wave  number,  wave  and  group  celerities,  and 
ray  angle  in  the  previous  column  are  calculated.  The  energy  is  calculated  as  a 
»veighted  average  of  energy  between  the  two  adjacent  grid  points  in  the  column 
and  the  direction  bins.  The  energy  density  is  corrected  by  a  factor  that  is  the 
ratio  of  the  5-deg  standard  angle  band  width  to  the  width  of  the  back-traced  band 
to  account  for  the  different  angle  increment  in  the  back-traced  ray.  The  shoaled 
and  refracted  wave  energy  in  column  2  is  then  calculated  from  the  conservation 
of  wave  action  along  a  ray  (Equation  9). 

In  a  strong  opposing  current  (e.g.,  ebb  currents  at  an  inlet),  waves  may  be 
blocked  by  the  current.  Blocking  occurs  if  there  is  no  solution  to  the  dispersion 
equation  (Equation  2).  Or,  to  state  it  another  way,  blocking  occurs  if  the  relative 
wave  group  celerity  is  smaller  than  the  magnitude  of  the  opposing  current,  so 
wave  energy  cannot  propagate  against  the  current.  In  deep  water,  blocking 
occurs  for  an  opposing  current  with  magnitude  greater  than  one-fourth  the 
deepwater  wave  celerity  without  current  (0.25  g  TJQ.tC),  where  is  the  absolute 
wave  period).  If  blocking  occurs,  the  wave  ener^  is  dissipated  through  breaking. 


Diffraction 

Diffraction  is  included  in  STWAVE  in  a  simple  manner  through  smoothing 
of  wave  energy.  The  model  smooths  energy  in  a  given  frequency  and  direction 
band  using  the  following  form 

£,(w„,e)  =  o.55£.(w^,e)  +  0.225  (£.^(w„,e)  +  £;.,(w„,e))  (lo) 


where  E  is  the  energy  density  in  a  given  frequency  and  direction  band,  and  the 
subscript  j  indicates  the  grid  row  index  (alongshore  position).  Equation  10  pro¬ 
vides  smoothing  of  strong  gradients  in  wave  height  that  occur  in  shelter  regions, 
but  provides  no  turning  of  the  waves.  This  formulation  is  grid-spacing  depen¬ 
dent,  which  is  a  serious  weakness.  Efforts  are  ongoing  to  implement  a  more 
rigorous  diffraction  representation. 


Source/sink  terms 

Surf-zone  wave  breaking.  The  wave-breaking  criterion  applied  in  the  first 
version  of  STWAVE  was  a  function  of  the  ratio  of  wave  height  to  water  depth 

H 

=  0.64  (11) 

d 
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where  H„,„  is  the  energy-based  zero-moment  wave  height.  At  an  inlet,  where 
waves  steepen  because  of  the  wave-current  interaction,  wave  breaking  is 
enhanced  because  of  the  increased  steepening.  Smith,  Resio,  and  Vincent 
(1997)  performed  laboratory  measurements  of  irregular  wave  breaking  on  ebb 
currents  and  found  that  a  breaking  relationship  in  the  form  of  the  Miche  criterion 
(1951)  was  simple,  robust,  and  accurate 

H  =  O.IZ  tanh  kd  (12) 

max 


(see  also  Battjes  1982  and  Battjes  and  Janssen  1978).  Equation  12  is  applied  in 
version  2.0  of  STWAVE  as  a  maximum  limit  on  the  zero-moment  wave  height. 
The  energy  in  the  spectrum  is  reduced  at  each  frequency  and  direction  in  propor¬ 
tion  to  the  amount  of  prebreaking  energy  in  frequency  and  direction  band.  Non¬ 
linear  transfers  of  energy  to  high  frequencies  that  occur  during  breaking  are  not 
represented  in  the  model. 

Wind  input.  Waves  grow  through  the  transfer  of  momentum  from  the  wind 
field  to  the  wave  field.  The  flux  of  energy,  into  the  wave  field  in  STWAVE 
is  given  by  (Resio  1988a) 

0.85C.  4  (13) 


where 

A  =  partitioning  coefficient  that  represents  percentage  of  total  atmosphere 
to  water  momentum  transfer  that  goes  directly  into  wave  field  (0.75) 


=  density  of  air 
p„,  =  density  of  water 
C„,  =  mean  wave  celerity 

M.  =  friction  velocity  (equal  to  product  of  wind  speed,  U,  and  square  root  of 
drag  coefficient,  Cd  =  .0012+.000025t/) 

In  deep  water,  STWAVE  provides  a  total  energy  growth  rate  that  is  consistent 
with  Hasselmann  et  al.  (1973). 

The  energy  gain  to  the  spectrum  is  calculated  by  multiplying  the  energy  flux 
by  the  equivalent  time  for  the  wave  to  travel  across  a  grid  cell 


At  = 


Ax 

PC^  cos 


(14) 
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where 


A?  =  equivalent  ttavel  time 
Ax  =  grid  spacing 


p  =  factor  equal  to  0.9  for  wind  seas 

Cg  =  average  group  celerity  of  spectrum 

0„  =  mean  wave  direction,  relative  to  grid 

Because  STWAVE  is  a  half-plane  model,  only  winds  blowing  toward  the  shore 
(+x  direction)  are  included.  Wave  damping  by  offshore  winds  and  growth  of 
offshore-traveling  waves  are  neglected. 


Wave-wave  interaction  and  white  capping.  As  energy  is  fed  into  the  waves 
from  the  wind,  it  is  redistributed  through  nonlinear  wave-wave  interaction. 
Energy  is  transferred  from  the  peak  of  the  spectrum  to  lower  frequencies 
(decreasing  the  peak  frequency  or  increasing  the  peak  period)  and  to  high 
frequencies  (where  it  is  dissipated). 


In  STWAVE,  the  frequency  of  the  spectral  peak  is  allowed  to  increase  with 
fetch  (or  equivalently  propagation  time  across  a  fetch).  The  equation  for  this  rate 
of  change  of  ^  is  given  by 


= 


<f^T  - 1  c 


/  \ 
u 


4/3 


At 


-3/7 


(15) 


where  the  i  and  i+l  subscripts  refer  to  the  grid  column  indices  within  STWAVE 
and  C  is  a  dimensionless  constant  (Resio  and  Perrie  1989).  The  energy  gained  by 
the  spectrum  is  distributed  within  frequencies  on  the  forward  face  of  the  spectrum 
(frequencies  lower  than  the  peak  frequency)  in  a  manner  that  retains  the  self¬ 
similar  shape  of  the  spectrum. 

Wave  energy  is  dissipated  (most  notably  in  an  actively  growing  wave  field) 
through  energy  transferred  to  high  frequencies  and  dissipated  through  wave 
breaking  (white  capping)  and  turbulent/viscous  effects.  There  is  a  dynamic  bal¬ 
ance  between  energy  entering  the  wave  field  because  of  wind  input  and  energy 
leaving  the  wave  field  because  of  nonlinear  fluxes  to  higher  frequencies  (Resio 
1987,  1988a).  The  energy  flux  to  high  frequencies  is  represented  in  STWAVE  as 


1/2  r'S  j 

eg  kp 
tanh^^'*  (k^  d) 


(16) 
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where 


/^  =  energy  flux 
e  =  coefficient  equal  to  30 
E  =  total  energy  in  spectrum 

kp  =  wave  number  associated  with  peak  of  spectrum  (Resio  1987) 

The  energy  loss  to  tiie  spectrum  is  calculated  by  multiplying  the  energy  flux  by 
the  equivalent  time  for  the  wave  to  travel  across  a  grid  cell  (At)  with  p  equal  to 
1.0  for  the  swell  portion  of  the  spectrum  and  0.9  for  the  sea  portion  of  the  spec¬ 
trum.  This  dissipation  is  only  applied  in  the  model  if  wind  input  is  included. 


Numerical  Discretization 

STWAVE  is  a  finite-difference  numerical  model,  formulated  on  a  Cartesian 
grid.  Grid  cells  are  square  (Ax  =  Ay).  Variable  grid  resolution  can  be  obtained 
by  nesting  model  runs.  This  is  accomplished  by  running  the  model  at  a  coarse 
resolution  and  saving  a  spectrum  at  a  nearshore  point  This  nearshore  spectrum 
can  then  be  used  as  a  boundary  condition  for  another  grid  of  finer  resolution.  A 
schematic  of  a  grid  is  shown  in  Figure  3.  STWAVE  operates  in  a  local  coordinate 
system,  with  the  x  axis  oriented  in  the  cross-shore  direction  (origin  offshore)  and 
tire  y  axis  oriented  alongshore,  forming  a  right-handed  coordinate  system.  The 
orientation  of  the  x  axis  (±87.5  deg)  defines  the  half  plane  that  is  represented  in 
the  model.  The  y  axis  is  typically  aligned  with  the  bottom  contours.  Wave  angles 
are  defined  in  a  mathematical  sense,  measured  counterclockwise  from  the  x  axis. 

Lateral  boundaries  in  the  model  can  be  specified  as  land  or  water  by  specify¬ 
ing  the  cell  depths  as  positive  (water)  or  negative  (land).  Note  that  specifying 
lanH  around  the  entire  grid  will  give  different  results  than  if  the  lateral  boundaries 
are  water.  Land  boundaries  reduce  wave  growth  near  the  boundary  because  they 
“block”  propagation  firom  landward  directions.  If  the  boundaries  are  specified  as 
water,  a  zero-gradient  type  of  boundary  is  applied  that  allows  energy,  consistent 
with  neighboring  cells,  to  propagate  into  or  out  of  the  domain  along  the  lateral 
boundary. 
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Figure  3.  Schematic  of  STWAVE  grid 
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Input  File  Description 


STWAVE  has  four  input  files.  These  files  specify  model  parameters, 
bathymetiy,  incident  wave  spectra,  and  current  fields.  The  current  field  file  is 
optional  and  required  only  if  the  wave-current  interaction  is  specified.  The  other 
three  files  are  always  required.  The  input  files  are  most  easily  generated  using 
the  SMS  interface  discussed  in  Chapter  5,  but  may  also  be  generated  using  an 
editor  or  other  programs.  File  formats  for  the  input  files  must  remain  the  same, 
regardless  of  how  the  files  are  generated.  Parameter  names  starting  with  “N”  or 
“I”  indicate  integer  parameter  values  are  required.  All  other  names  indicate  real 
parameter  values.  A  schematic  of  the  STWAVE  input  and  output  files  is  given 
in  Figure  4. 


Figure  4.  STWAVE  input  and  output  files 
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Model  Parameters 


The  model-parameter  file  specifies  options  for  running  STWAVE  and  special 
output  points.  The  first  line  of  the  model  parameter  file,  read  using  free  format 
(parameters  must  be  delimited  by  spaces,  commas,  or  new  lines),  includes  (in  the 
following  order): 

IPRP  =  Switch  for  including  propagation  only  (IPRP  =  1)  or  both  propagation 
and  source  terms  (IPRP  =  0).  Surf-zone  wave  breaking  is  included  for  either 
option,  but  wind-wave  generation,  wave-wave  interactions,  and  white  capping 
are  included  only  for  IPRP  =  0.  In  applications  where  the  wave  propagation 
distances  are  short  (on  the  order  of  kilometers),  propagation  only  is  usually 
sufficient.  For  longer  wave  propagation  distances  or  cases  where  only  locally 
generated  waves  exist,  the  source  terms  should  be  included. 

ICUR  =  Switch  for  including  (ICUR  ==  1)  or  excluding  (ICUR  =  0)  wave- 
current  interaction.  A  current  field  must  be  specified  (see  Current  Fields  below) 
if  ICUR  =  1 .  Currents  can  increase  (opposing  current)  or  decrease  (following 
current)  wave  height  as  well  as  change  wave  direction  because  of  wave-current 
interaction.  These  interactions  are  strongest  for  shore  wave  periods  and  large 
current  magnitudes.  Smith,  Militello,  and  Smith  (1998)  showed  that  current 
magnitudes  less  than  about  1  m/sec  had  minor  impacts  on  wave  transformation  at 
an  East  Coast  inlet  in  Florida  (peak  wave  periods  in  the  range  of  5-1 5  sec). 

NSELCT  =  Number  of  special  output  points.  The  wave  spectra  will  be  saved 
only  for  these  output  points.  Also,  a  file  that  summarizes  wave  height,  period, 
and  direction  for  these  points  will  be  written.  Wave  height,  period,  and  direction 
are  saved  by  the  model  at  all  grid  points  (see  Chapter  4),  but  this  summaiy  file 
makes  it  easy  to  spot  check  results  at  critical  locations  or  to  isolate  only  the 
locations  of  interest. 

The  final  information  in  the  model-parameter  file  is  specification  of  the 
special  output  points.  These  points  are  I,  J  {x  grid  cell  index,  y  grid  cell  index) 
pairings  of  these  special  output  points.  There  must  be  NSELCT  pairs  of 
numbers.  The  values  are  again  read  in  free  format  and  may  occupy  as  many 
lines  as  needed. 

The  default  name  of  the  model  parameter  files  is  options .  std.  In  SMS, 
the  file  is  specified  with  a  project  name  and  the  extension  .  std.  A  sample 
model  parameter  file  is  give  in  Appendix  A.  The  sample  file  specifies  that  the 
model  run  will  consider  only  propagation  with  no  source  terms  (IPRP  =  1),  and 
wave-current  interaction  is  included  (ICUR  =  1).  Three  output  points  have  been 
selected  (at  I,  J  coordinates  (15,21),  (15,65),  and  (28,49)). 
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Bathymetry 

The  bathymetry  file  gives  the  STWAVE  grid  size  and  spacing  and  defines  the 
grid  bathymetiy.  The  first  line  of  the  file,  read  using  free  format,  includes  (in  the 
following  order): 

NI  =  Number  of  cross-shore  (x)  grid  cells  or  columns.  The  value  of  NI, 
together  with  the  grid  spacing  (DXINC),  determines  the  cross-shore  extent  of  the 
modeling  domain  and  the  location  of  the  offshore  grid  boundaiy.  The  offshore 
boundary  is  typically  placed  in  a  region  where  the  bathymetry  contours  are 
relatively  straight  and  parallel  and  where  incident  wave  information  is  available 
from  a  gauge  or  hindcast. 

NJ  =  Number  of  alongshore  (y)  grid  cells  or  rows.  The  value  of  NJ,  together 
with  the  grid  spacing,  determines  the  alongshore  extent  of  the  modeling  domain 
and  the  location  of  the  lateral  grid  boundaries.  The  lateral  extent  of  the  grid  is  a 
balance  between  keeping  the  lateral  boundaries  as  far  away  as  possible  from  the 
region  of  engineering  interest  and  keeping  the  grid  domain  small  for  computa¬ 
tional  efficiency.  The  grid  should  extend  sufficiently  far  alongshore  to  include 
the  bathymetric  features  that  influence  the  engineering  project  (e.g.,  if  a  project 
is  centered  on  a  coastal  inlet  entrance,  the  entire  ebb  shoal  and  any  offshore  fea¬ 
tures  that  influence  refraction  should  be  included  in  the  model  domain). 

DXINC  =  Grid  spacing  in  meters,  which  is  the  same  in  the  x  and  y  directions. 
The  grid  spacing  determines  how  finely  bathymetric  features  and  current  fields 
are  resolved  in  the  model  and  the  spatial  resolution  of  the  output  fields  of  wave 
height,  period,  and  direction.  Grid  spacing  should  be  sufficient  to  define  bathy¬ 
metric  features  and  gradients  (shoals,  canyons,  and  channels)  that  are  important 
to  the  engineering  project.  Typical  resolutions  are  tens  to  hundreds  of  meters. 
For  a  fixed  grid  domain,  increasing  the  grid  resolution  increases  computational 
time. 

The  remainder  of  the  file  gives  the  water  depth  (relative  to  some  datum,  such  as 
mean  sea  level  or  mean  lower  low  water)  for  each  grid  cell  in  an  STWAVE 
application,  in  meters.  Water  depths  are  defined  as  positive  numbers  and  land 
elevations  are  negative  numbers.  The  depths  are  read  with  a  free  format.  The 
file  begins  with  the  depth  at  cell  (1 ,  NJ)  and  reads  in  the  cross-shore  direction 
(I  =  1  to  NI).  The  read  is  repeated  for  J  =  NJ-1  (I  =  1  to  NI),  and  progresses  to 
J  =  1 .  The  default  name  for  the  bathymetry  file  is  depth .  in.  In  SMS,  the  file 
is  specified  with  a  project  name  and  the  extension  .  dep.  A  sample  bathymetry 
file  is  given  in  Appendix  B.  The  sample  grid  is  defined  as  53  cells  in  the  cross¬ 
shore  direction  and  1 12  cells  in  the  alongshore  direction.  The  grid  spacing  is 
50  m. 
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Incident  Wave  Spectra 

Incident  two-dimensional  wave  spectra  are  specified  as  energy  density  as  a 
function  of  frequency  and  direction.  A  single  input  spectrum  is  applied  along 
the  entire  offshore  boundary  of  the  STWAVE  grid  (the  spectrum  is  set  to  zero 
for  any  land  points  along  the  boundaiy).  Thus,  it  is  good  practice  to  establish  the 
offshore  boundaiy  of  the  bathymetiy  grid  along  a  depth  contour  where  the  wave 
spectrum  is  fairly  homogeneous  (no  large  shoals  or  canyons  offshore  of  the 
boundary).  The  first  line  of  the  spectra  file  describes  the  number  of  frequency 
and  direction  bins  in  the  spectra  (read  in  free  format): 

NF  =  Number  of  frequency  bins  in  the  spectra.  The  number  of  frequency 
bins  determines  how  finely  resolved  are  the  calculated  spectra.  A  large  number 
of  frequency  bins  increases  the  computation  time,  and  a  small  number  of  bins 
reduces  model  resolution.  Typically,  20-30  bins  are  used. 

NA  =  Number  of  direction  bins  in  the  spectra.  This  value  must  be  set  to  35, 
which  gives  5-deg  resolution  in  direction. 

The  next  lines  of  the  model  parameter  file  specify  the  frequencies  for  model 
spectra  (used  for  the  input  spectra,  internal  computations,  and  output  spectra), 
starting  from  the  lowest  frequency.  There  must  be  NF  frequencies  specified. 

The  frequencies  are  again  read  in  free  format  and  may  occupy  as  many  lines  as 
needed.  These  frequencies  should  span  the  frequency  range  where  significant 
wave  energy  is  contained  in  the  spectrum.  This  can  be  estimated  by  inspecting 
the  input  spectrum  or  estimating  the  peak  period  expected  using  the  wave  growth 
curves  in  the  Coastal  Engineering  Manual  (in  preparation).  A  rule  of  thumb  is 
that  the  spectral  peak  should  fall  at  about  the  lower  one-third  of  the  frequency 
range  (e.g.,  if  the  peak  frequency  is  0.1  Hz,  the  range  may  be  0.01  to  0.3  Hz). 
Wave  frequencies  higher  (or  periods  shorter)  than  the  highest  frequency  bin  or 
lower  than  the  lowest  frequency  bin  will  not  be  resolved  by  the  model.  Typi¬ 
cally,  frequency  increments  are  on  the  order  of  0.01  Hz,  but  the  increment  need 
not  be  constant.  West  Coast  applications  will  tend  to  require  finer  resolution 
focused  at  lower  frequencies  (because  of  long  wave  periods),  and  Gulf  Coast  or 
Great  Lakes  applications  will  tend  to  require  coarser  resolution  covering  a 
broader  range  of  frequencies  (because  of  shorter  wave  periods). 

Following  specification  of  the  frequency  and  direction  bins  is  a  header  line 
containing  a  spectrum  identifier,  wind  information,  peak  frequency,  and  water- 
elevation  correction: 

IDD  =  Integer  event  identifier,  such  as  a  date  or  test  number.  The  integer  is 
limited  to  nine  digits. 

U  =  Wind  speed  in  meters/second.  This  wind  speed  is  considered  constant 
over  the  entire  grid.  Winds  that  are  blowing  offshore  (-x  direction)  are  neglected 
by  the  model,  because  the  model  only  treats  generation  in  a  half  plane.  For  wind 
directions  greater  than  60  deg  relative  the  x  axis,  STWAVE  underpredicts  wave 
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generation  (because  of  the  half-plane  model  coverage).  To  model  large  ranges  of 
wind  direction,  multiple  STWAVE  grid  orientations  may  be  required. 

UDIR  =  wind  direction  relative  to  the  STWAVE  coordinate  system  in 
degrees.  Wind  direction  is  measured  counterclockwise  from  the  x  axis  (same 
convention  used  for  wave  direction,  see  Figure  3). 

FM  =  Peak  spectral  frequency  in  Hz.  If  a  spectrum  is  specified  with  no 
energy  (i.e.,  waves  will  be  locally  generated  in  the  modeling  domain),  FM  should 
be  set  to  the  highest  frequency  specified  in  the  model  parameter  file. 

DADD  =  Water-elevation  correction  (tide  or  storm  surge)  in  meters,  relative 
to  the  bathymetry  datum.  Users  should  take  special  care  to  ensure  a  consistent 
datum  is  used  to  specify  the  bathymetry  and  the  water-elevation  correction.  The 
water-elevation  correction  is  constant  over  the,  entire  grid. 

The  header  line  is  read  in  free  format  and  is  followed  by  the  energy  densities 
in  the  units  meters  squared/hertz/radian.  The  spectrum  is  read  starting  with  the 
lowest  frequency  and  reading  all  the  wave  directions  (from  -85  deg  to  85  deg), 
then  reading  energy  densities  for  all  directions  for  the  next  lowest  frequency,  etc. 
The  total  number  of  energy  densities  read  equals  the  number  frequencies  times 
the  number  of  directions. 

Multiple  wave  model  runs  can  be  executed  by  inserting  multiple  spectra,  with 
identifying  header  lines,  in  the  input  file  (the  number  of  frequency  and  direction 
bins  and  the  frequencies  are  not  repeated  for  each  spectrum).  The  model  will 
loop  through  the  execution  until  it  runs  out  of  input  spectral  information  (model 
parameters  and  bathymetry  will  stay  constant).  If  wave-current  interaction  is 
specified,  the  number  of  current  fields  in  the  current  input  files  must  equal  the 
number  of  input  spectra.  The  default  name  for  the  incident  spectra  file  is 
spec .  in.  In  SMS,  the  file  is  specified  with  a  project  name  and  the  extension 
.  eng  (for  energy  density).  A  sample  incident  spectrum  file  is  given  in 
Appendix  C.  The  sample  spectrum  file  has  30  frequencies  and  35  directions. 

The  spectral  frequencies  are  from  0.031  to  0.258  Hz,  incrementing  by  0.0078  Hz. 
The  run  identifier  is  96031202  (12  March  1996  at  0200  GMT).  The  wind  speed 
and  wind  direction  are  zero  because  local  wind  generation  is  not  specified  in  the 
run  (IPRP  =  1 ,  in  the  model  parameter  file).  The  peak  frequency  of  the  spectrum 
is  0.078  Hz  (period  of  13  sec),  and  the  tide/surge  elevation  is  0.34  m  relative  to 
mean  sea  level  (which  corresponds  to  the  mean  sea  level  datum  used  to  develop 
the  bathymetry  file). 

The  incident  wave  spectrum  can  be  specified  using  a  measured  two- 
dimensional  spectrum,  by  generating  a  spectrum  using  a  parametric  shape,  or 
from  a  spectrum  calculated  by  a  global-  or  regional-scale  wave  model.  Rarely  do 
measured  spectra  have  the  directional  resolution  required  by  the  model;  thus, 
measured  spectra  require  interpolation  to  finer  resolution.  A  parametric  spectral 
shape  together  with  a  directional  spreading  function  can  also  be  applied  to 
specify  an  incident  spectrum,  knowing  the  wave  height,  period,  and  direction. 
SMS  includes  the  capability  to  generate  incident  spectra  using  a  TMA 
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one-dimensional  shallow-water  spectral  shape  (named  for  the  three  data  sets 
used  to  develop  the  spectrum:  TEXEL  storm,  MARSEN,  and  ARSLOE)  (Bouws 
et  al.  1985)  and  a  cos""0  directional  distribution  (see  Chapter  5).  To  generate  a 
TMA  spectrum,  the  following  parameters  must  be  specified:  peak  wave  period 
(Tp),  wave  height,  water  depth,  and  a  spectral  peakedness  parameter  (y).  The 
peakedness  parameter  controls  the  width  of  the  frequency  spectrum  (small 
numbers  give  broad  peaks  and  large  numbers  give  narrow  peaks).  The  direc¬ 
tional  distribution  of  the  spectrum  is  specified  with  a  mean  direction  and  a  direc¬ 
tional  spreading  coefficient  («w).  The  energy  in  the  frequency  spectrum  is 
spread  proportional  to  cos"”(6-6J,  where  0  is  direction  of  the  spectral  compo¬ 
nent  and  0„  is  the  mean  wave  direction.  Guidance  for  selecting  y  and  nn  is 
given  in  Table  1 Spectra  with  large  peak  periods  are  generally  narrow  in  both 
frequency  and  direction  (swell).  For  small  peak  periods,  spectra  are  typically 
broad  in  both  frequency  and  direction  (sea). 


Table  1 

Approximate  Spectral  Peakedness  and  Directional  Spreading 
Parameters 

r,,  Sec 

Y 

nn 

<10 

3.3 

4 

11 

4 

8 

12 

4 

10 

13 

5 

12 

14 

5 

16 

15 

6 

18 

16 

6 

20 

17 

7 

22 

18 

7 

26 

19 

8 

28 

20 

8 

30 

Current  Fields 

A  current  field  input  file  is  required  only  if  ICUR  =  1  in  the  model  parameters 
(it  is  ignored  if  ICUR  =  0).  The  current  vector  is  specified  as  x  and  y  compo¬ 
nents  of  the  current  at  each  grid  cell,  u  and  v,  in  the  units  meters/second.  The  x 
and  y  directions  correspond  to  the  x  and  y  axes  of  the  grid  in  the  local  STWAVE 
coordinates  (Figure  3).  The  first  line  of  the  current  input  file  is  the  Nl,  NJ,  and 


'  Personal  Communication,  1998,  Edward  F.  Thompson,  U.S.  Army  Engineer  Waterways 
Experiment  Station,  Vicksburg,  MS. 
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DXINC  repeated  from  the  bathymetiy  file.  STWAVE  checks  that  these  values 
are  consistent  with  the  bathymetiy  file,  and  if  they  are  not,  model  execution  is 
terminated  with  the  error  “current  field  does  not  match  depth  grid  size.”  Next, 
the  header  for  the  current  field  is  given,  which  is  an  event  identifier.  This  header 
is  not  used  by  the  model,  but  is  helpful  to  ensure  that  the  current  and  spectral 
files  coincide.  Similar  to  the  bathymetry  file,  the  current  file  begins  with  the  u/v 
pair  at  cell  (1,  NJ)  and  reads  in  the  cross-shore  direction  (I  =  1  to  NI).  Then  the 
read  is  repeated  for  J  =  NJ-1  (1  =  1  to  NI),  and  progresses  to  J  =  1 .  The  current 
components  are  read  in  free  format.  The  default  name  for  the  current  field  file  is 
current .  in.  In  SMS,  the  file  is  specified  with  a  project  name  and  the  exten¬ 
sion  .  cur.  A  sample  current  field  file  is  given  in  Appendix  D.  In  the  sample 
file,  the  first  line  repeats  the  grid  parameters  from  the  bathymetry  file.  The  grid 
is  defined  as  53  cells  in  the  cross-shore  direction  and  112  cells  in  the  alongshore 
direction  with  a  grid  spacing  of  50  m.  The  event  identifier  for  the  current  field  is 
96031202,  which  matches  the  event  identifier  in  the  incident  spectrum  file.  As 
with  the  spectral  file,  multiple  wave  model  runs  can  be  executed  by  inserting 
current  fields  (with  their  event-identifier  headers)  in  the  input  file  (as  well  as  the 
same  number  of  input  spectra  in  the  spec.in  file).  The  model  will  loop  through 
the  execution  until  it  runs  out  of  input  spectral  and  current  information. 
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4  Output  File  Description 


STWAVE  has  three  output  files.  These  files  contain  wave  spectra  at  the 
selected  output  points;  wave  height,  period,  and  direction  at  the  selected  output 
points;  and  fields  of  wave  height,  period,  and  direction  over  the  entire  STWAVE 
modeling  domain.  If  multiple  incident  wave  spectra  are  specified  in  the  input 
file  (and  current  fields,  if  wave-current  interaction  is  specified),  then  output  from 
each  of  these  model  runs  is  appended  to  the  output  files. 


Selected  Wave  Spectra 

Wave  spectra  at  selected  grid  cells,  specified  in  the  model  parameters  file 
(see  Chapter  3,  Input  File  Description),  are  saved  in  a  spectral  output  file. 

Similar  to  the  input  spectra  file,  the  first  line  of  the  output  file  describes  the 
number  of  frequencies  and  number  of  directions  (NF  and  NA),  and  the  following 
lines  give  the  frequencies.  Each  spectrum  in  the  file  is  proceeded  by  the 
following  header  information: 

IDD  =  Integer  event  identifier,  specified  in  the  incident  wave  spectrum  file. 

I  =  Integer  X  grid  cell  index  of  the  special  output  point,  specified  in  the  model 
parameter  file. 

J  =  Integer^  grid  cell  index  of  the  special  output  point,  specified  in  the  model 
parameter  file. 

N  =  Integer  index  of  the  special  output  point  (values  are  1  to  NSELCT,  based 
on  the  order  in  which  the  points  are  specified  in  the  model  parameter  file). 

The  header  format  is  (i  1 0,3i5).  Following  the  header  line,  the  output  spectrum  is 
written  in  the  same  order  and  units  in  which  the  incident  spectrum  was  read. 
Energy  densities  are  given  in  the  units  meters  squared/hertz/radian.  The 
spectrum  is  written  starting  with  the  lowest  frequency  and  writing  energy  density 
for  all  wave  directions  (from  -85  deg  to  85  deg),  followed  by  the  energy  densities 
for  all  directions  for  the  next  lowest  frequency,  etc.  The  format  of  the  energy 
densities  is  (17f6.3). 
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The  spectra  will  be  written  to  the  file  in  an  order  based  on  the  I  cell  index, 
starting  from  the  smallest  index  and  proceeding  to  the  largest  index.  This  is 
because  spectra  are  not  saved  over  the  entire  grid,  so  the  output  spectra  are 
written  to  the  file  as  the  grid  column  is  solved  and  then  are  overwritten  in 
memory.  Thus,  the  order  may  not  coincide  with  N  =  1,  NSELCT.  The  default 
file  name  for  the  output  spectra  file  is  spec .  out.  In  SMS,  the  file  is  specified 
with  a  project  name  and  the  extension  .obs. 


Selected  Wave  Parameters 

The  wave  parameters  H„,^  Tp,  and  mean  wave  direction,  0„„  at  selected  grid 
cells,  specified  in  the  model  parameters  file  (see  Chapter  3,  Input  File 
Description),  are  saved  in  a  selected  parameter  output  file.  The  selected  grid 
cells  are  the  same  locations  where  spectra  were  saved.  For  each  output  point,  the 
following  information  is  written: 

IDD  =  Integer  event  identifier,  specified  in  the  incident  wave  spectrum  file. 

I  =  Integer  X  grid  cell  index  of  the  special  output  point,  specified  in  the  model 
parameter  file. 

J  =  Integer  y  grid  cell  index  of  the  special  output  point,  specified  in  the  model 
parameter  file. 

=  Zero-moment  wave  height  at  cell  (I,  J),  in  meters. 

Tp  =  Peak  wave  period  at  cell  (I,  J),  in  seconds. 

0„,  =  Mean  wave  direction  at  cell  (I,  J),  in  degrees  relative  to  the  STWAVE 
grid  (see  Figure  3). 

The  format  of  each  line  is  (il0,2i3,f6.2,f6.1,f6.0).  Information  for  each  selected 
grid  cell  is  written  on  a  separate  line  in  the  file.  The  default  file  name  is 
selhts .  out.  This  file  is  not  used  in  SMS  because  wave  height,  period,  and 
direction  for  all  grid  cells  are  displayed  graphically  on  the  screen  over  the  entire 
domain  (values  at  any  cell  in  the  domain  are  displayed  by  clinking  on  the  cell 
with  the  mouse).  Examples  of  two  selected  wave  parameter  files  are  given  in 
Chapter  6. 


Wave  Parameter  Fields 

The  wave  parameters  H„,„  Tp,  and  are  also  saved  for  all  grid  cells  into  a 
wave  field  file.  Written  first  in  the  file  are  the  grid  dimensions  (NI  and  NJ)  and 
grid  spacing  (DXINC)  repeated  from  the  bathymetry  input  file.  Next  in  the  file 
is  the  event  identifier  (IDD  from  the  incident  wave  spectrum  file)  in  the  format 
(i  1 0).  After  this  header  line,  all  the  wave  heights  are  written,  beginning  at  cell 
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(1,  NJ)  and  proceeding  in  the  cross-shore  direction  (1=1,  NI).  The  write  is 
repeated  for  each  grid  row  from  J  =  NJ- 1  to  J  =  1 .  The  format  is  ( 1 6f5.2). 
Following  the  wave  heights,  the  wave  periods  are  written  in  the  same  order  with 
the  format  (16f5.1).  Then  the  wave  directions  are  written  (again  in  the  same 
order)  with  the  format  (16i5).  If  multiple  input  spectra  are  specified,  the  event 
identifier  and  the  wave  heights,  periods,  and  directions  for  each  condition  are 
appended  to  the  end  of  the  file. 

The  default  file  name  for  the  wave  field  file  is  wavfid.  In  SMS,  the  file  is 
specified  with  a  project  name  and  the  extension  .wav. 
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5  SMS  User  Interface 


Information  contained  in  the  STWAVE  input  files  (Chapter  3)  and  output 
files  (Chapter  4)  is  conceptually  easy  to  understand,  but  practically,  these  files 
contain  large  amounts  of  data  that  can  be  difficult  to  generate  and  interpret.  To 
assist  STWAVE  users  in  generating  input  files  and  visualizing  output  files,  a 
user  interface  has  been  built  for  STWAVE  within  the  Surface-water  Modeling 
System  (Brigham  Young  University  Engineering  Computer  Graphics  Laboratory 
(ECGL)  1997).  The  SMS  interface  supports  grid  generation,  interpolation  of 
current  fields,  generation  of  input  spectra,  visualization  of  wave  heights,  periods, 
and  directions,  and  visualization  of  output  spectra.  STWAVE  will  be  released 
SMS  version  6.1  in  early  1999.  Corps  of  Engineers  SMS  users  are  supported 
through  the  Coastal  and  Hydraulics  Laboratory 

(http://hlnet.wes.army.mil/software/sms/),  and  non-Corps  users  can  get 
information  from  ECGL  (http://www.ecgl.byu.edu).  SMS  is  not  required  to  run 
STWAVE,  but  it  is  a  powerful  tool  to  efficiently  develop  and  visualize  input  and 
output  files. 

This  chapter  provides  a  brief  tutorial  on  the  application  of  STWAVE  within 
SMS.  It  does  not  cover  all  the  SMS  features,  but  describes  a  typical  application 
of  the  model.  For  coastal  inlets,  it  is  common  to  run  the  finite-element 
circulation  model  ADCIRC  (Luettich,  Westerink,  and  Scheffner  1992)  to 
calculate  water-surface  elevations  and  current  fields,  as  well  as  to  run  STWAVE. 
The  example  application  starts  with  an  ADCIRC  finite-element  bathymetry  grid 
(file  ponce.grd)  and  a  current  field  generated  by  ADCIRC  (file  ponce.64). 
These  gridded  data  are  converted  to  “scatter”  data  (random  data  points  that  are 
not  associated  with  a  grid,  but  with  x  andy  locations).  From  these  scatter  data,  a 
Cartesian  bathymetry  grid  and  current  field  are  generated.  Next,  other  required 
model  input  (input  spectrum  and  model  parameters)  is  entered  and  the  model  is 
executed.  Finally,  postprocessing  of  the  data  is  described,  including  visualiza¬ 
tion  of  the  bathymetry  and  wave  heights,  directions,  heights,  and  spectra. 
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Converting  ADCIRC  to  Scatter 

Reading  in  the  ADCIRC  files 

1.  Go  to  ihQ  Mesh  Module  and  choose,  Data\Switch  Current  Model...  Click 
on  ADCIRC  and  push  OK. 

2.  Go  to  ADCIRC\Open  Simulation...  Toggle  the  Control  (*.  1 5)  button  off 
and  make  sure  the  Open  Grid  File  button  is  on.  Click  the  button  under 
Open  Grid  File  and  choose  ponce .  grd.  Push  OK  to  read  in  the  mesh 
file. 

3.  Go  to  Data\Data  Browser  and  click  the  Import...  button  at  the  bottom 
left.  Change  the  F//e  7>/7e  to  ADCIRC  and  push  OK.  Select  ponce.  64 
and  push  OK. 

4.  After  the  data  is  read  in,  push  the  Convert  button  (still  in  the  Data 
Browser)  in  the  Vector  Data  Set  section.  Choose  Vx  and  Vy  to  split  the 
vector  data  into  its  components.  Rename  the  Vx  and  Vy  data  sets  to 
“Current_x”  and  “Current__y”  and  push  Convert.  When  text  appears 
telling  you  that  the  conversion  is  finished,  push  Done. 

5 .  Push  Done  again  to  exit  the  Data  Set  Browser. 


Converting  to  scatter 

1 .  Go  to  Data\Mesh  ->  Scatterpoint,  leave  the  name  as  “Scatter,”  and  push 
OK. 

2.  Switch  to  the  Scatter  Module  and  zoom  in  on  the  dense  part  as  shown  in 
Figure  5.  Choose  Edit\Select  With  Poly  and  click  on  the  screen  to  select 
enough  scatter  points  to  surround  the  area  where  the  grid  will  be  created. 
Double  click  on  the  screen  to  close  the  polygon  and  select  the  points. 

3.  Go  to  Interpolation\Create  Subset.  Rename  the  scatter  set  to  “Grid”  and 
push  OK. 

4.  Push  the  Frame  button  to  show  both  scatter  sets. 

5.  Select  the  Select  Scatter  Sets  tool  and  click  in  the  box  on  the  screen 
belonging  to  the  larger  scatter  set  called  “Scatter”  and  push  the  DELETE 
button.  This  deletes  the  unnecessary  data. 
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Figure  5.  Zoomed-in  view  of  scatter  data  set 


Creating  the  Cartesian  Grid 

Creating  the  Cartesian  grid  frame 

1 .  Switch  to  the  Map  Module.  Go  to  Feature  Objects\Coverages. . .  and 
switch  the  Coverage  Type  to  Cartesian  Grid.  Push  OK. 

2.  Go  to  Feature  Objects  \  Grid  Frame  and  push  the  New  Grid  button. 

3.  Drag  and  resize  the  grid  frame  by  dragging  the  comers  or  edges  until  the 
grid  frame  fits  over  the  desired  area.  Dragging  a  comer  resizes  the  frame. 
Dragging  an  edge  moves  the  entire  frame. 

4.  Rotate  the  grid  by  dragging  the  circle  at  the  bottom  right  comer  of  the 
frame.  Rotate  the  frame  until  the  circle  is  at  the  left  of  the  screen.  The 
origin  of  the  grid  frame  is  the  frame  comer  at  the  top  of  the  screen.  Note: 
It  may  be  necessary  to  make  the  grid  frame  smaller  and  then  push  OK  to 
exit  the  dialog.  You  can  zoom  in  on  the  desired  area  and  then  position 
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the  grid  frame  by  reopening  the  grid  frame  dialog  with  Feature 
Objects\Grid  Frame. 


5.  Push  OK  when  you  are  done  positioning  the  frame  to  appear  as  Figure  6. 


Figure  6.  Grid  frame  positioned  over  scatter  data  set 


Creating  the  land  and  ocean  polygons 

1 .  Go  back  to  the  Mesh  Module  and  choose  Data\Material-> Feature  and 
use  the  active  coverage.  This  creates  map  module  polygons  around  the 
boundaries  of  the  mesh. 

2.  Go  to  Display\Display  Options...  and  click  All  off  and  OK  to  turn  off  the 
elements  and  nodes. 

3.  Switch  back  to  the  Map  Module. 

4.  Select  the  Create  Arcs  tool  and  create  arcs  surrounding  the  left  part  of  the 
grid  frame  (outside  of  the  polygons  created  in  Step  1 ).  Connect  the  arcs 
with  vertices  of  the  ocean  polygon  as  shown  in  Figure  7. 
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Figure  7.  Arc  created  around  border  of  grid  frame 

5.  Go  to  Feature  Objects\Build Polygons  and  push  OK  at  the  prompt. 

6.  Double  click  in  the  new  polygon  that  you  created  (on  the  left  side)  and 
select  Land  and  push  OK.  The  default  for  the  polygons  created  earlier  is 
Ocean. 
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Mapping  to  the  grid 

1 .  Go  to  Feature  Objects\Map->2D  Grid. 

2.  Set  the  Cell  Options  to  Number  of  Columns  and  set  the  value  to  50. 

3.  Toggle  Create  Functions  on. 

4.  Click  the  button  that  is  to  the  right  of  Depth.  This  opens  the  Scatter 
Options  dialog. 

5.  Select  the  function  (“elevation”)  that  will  be  used  to  interpolate  the  depth 
and  push  OK. 

6.  Click  the  button  that  is  to  the  right  of  Current  x.  Select  the  function 
“Current  x”.  This  will  be  used  to  interpolate  the  X  component  of  current. 

7.  Change  Time  Steps  to  Multiple  and  change  Time  Step  Interpolation  to 
Interpolate. 

8.  Enter  the  26 1 000  in  Beg  time,  900  in  Step  size  and  3  in  Number  steps. 
This  tells  SMS  to  create  current  data  at  time  261000, 261900,  262800, 
and  263700  sec  (four  sets  of  values  around  three  steps  of  15  min  each). 

9.  Select  OK  to  exit  the  Scatter  Options  dialog. 

10.  Click  the  button  that  is  to  the  right  of  Current y.  Select  the  function 
“Current_y”.  This  will  be  used  to  interpolate  the  Y  component  of  current. 
You  do  not  need  to  repeat  the  time  step  definition,  since  both  components 
share  the  same  parameters. 

11.  Select  OK  to  exit  the  Scatter  Options  dialog. 

12.  Push  OK  to  create  the  Cartesian  grid. 

13.  Switch  to  the  Cartesian  Grid  module  to  edit  and  look  at  the  completed 
grid. 

Note  on  interpolation  (When  interpolating  you  can  specify  a  single  time  step  or 
multiple  steps.  Single  times  can  come  from  any  time  in  the  data  set  or 
interpolated  between  two  of  these  times.  For  multiple  steps,  you  can  specify  to 
match  the  steps  from  the  data  set  for  a  specified  number  of  steps,  or  you  can 
specify  a  new  time  step  size  and  the  number  of  steps.  SMS  will  generate  one  set 
of  values  at  the  beginning  time  specified,  and  one  at  the  end  of  each  time  step. 
Therefore,  there  will  be  “n+l”  sets  of  currents  if  you  specify  “n”  time  steps.) 
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Editing  the  Grid  and  Running  STWAVE 

Generating  spectral  energy  distribution 

1 .  In  the  Cartesian  Grid  module,  choose  STWA  VE\Spectral  Energy. 

2.  Push  the  Generate  Spectrum  button.  Enter  40  for  Num.  Frequencies, 

0.0 1  for  the  Frequency  Step,  and  0.03  for  the  Starting  Frequency. 

3.  Push  the  Generate  button  and  then  OK  to  exit  the  Spectral  Energy  dialog. 

Model  control 

1.  Se.\tciSTWAVE\ModelControl. 

2.  Change  the  Wind  to  Propagation  Only  and  push  OK  to  exit  the  dialog. 

Selecting  monitoring  stations 

1 .  Select  the  Select  Cell  tool. 

2.  Select  the  cell  with  the  IJ  location  of  I  =  34  and  J  =  30.  The  location  can 
be  seen  in  the  bottom  right  in  the  Edit  Window  when  a  cell  is  selected. 
(Click  on  the  screen  to  select  a  cell.) 

3.  While  the  cell  is  still  selected,  choose  Ce// Turn 

the  Monitoring  Station  toggle  on  and  push  OK. 

4.  Repeat  Steps  2  and  3  and  change  the  cells  at  I  =  1 7,  J  =  1 8  and  I  =  25,  J  = 
35  to  be  Monitoring  cells. 

Saving  the  simulation 

To  save  the  simulation,  select  STWA  VE\Save  Simulation  and  enter  ponce .  sim. 


Running  STWAVE 

1 .  To  run  STWA  VE,  select  STWA  VE\Run  STWA  VE. 

2.  Click  the  folder  button  and  path  to  and  select  the  STWA  VE  executable. 
Push  OK. 

3.  A  window  will  appear  and  stay  up  while  runs.  Vlhen  STWAVE 

is  finished,  the  window  will  disappear. 
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Postprocessing 

Visualizing  the  STWAVE  solution 

1 .  Open  the  solution  files  by  selecting  STWA  VE\Open  Simulation. 

2.  Select  ponce .  s  im  and  push  OK. 

3.  Select  Display\Display  Options.  Turn  the  Contours  and  Vectors  toggles 
on. 

4.  Push  the  Contour  Options  button  and  select  Color  Fill  Between  Contours 
as  the  Contour  Method.  Push  OK. 

5 .  Push  the  Vector  Options  button  and  change  the  Shaft  Length  to  Define 
min  and  max. 

6.  Set  the  Min  length  to  25  and  the  Max  length  to  50.  Push  OK. 

7.  Push  OK  to  exit  the  Display  Options  dialog. 

Visualizing  bathymetry 

1 .  Go  to  Data\Data  Browser  and  select  the  scalar  function  named  “depth.” 
Push  Done  to  exit  the  dialog. 

2.  The  contours  show  the  water  depth  in  the  ocean. 

Visualizing  the  direction  field 

1 .  Go  to  Data\Data  Browser  and  select  the  vector  function  named  “wave.” 
Push  Done  to  exit  the  dialog. 

2.  The  vector  arrows  show  the  wave  direction  field. 

Visualizing  the  wave  height 

1 .  Go  to  Data\Data  Browser  and  select  the  scalar  function  named  “Height.” 
Push  Done  to  exit  the  dialog. 

2.  The  contours  show  the  wave  height. 
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Visualizing  the  observation  (monitoring  cells)  spectra 

1 .  Go  to  SUVA  VE\Spectral  Energy. 

2.  Push  the  Contour  Options  button,  select  Color  Fill  Between  Contours, 
and  push  OK. 

3.  Push  the  Data  Browser  button  and  select  the  scalar  data  set  named 
“Node_34_30.”  Push  Done. 

4.  The  contours  show  the  observation  spectra  for  the  node  at  I  =  34,  J  =  30. 

5.  Push  OK  to  exit  the  Spectral  Energy  dialog. 
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6  Example  Applications 


In  this  chapter,  two  example  applications  of  STWAVE  are  given.  The 
purpose  of  these  examples  is  to  demonstrate  how  the  model  can  be  applied  and 
provide  sample  problems  that  a  new  user  can  follow  to  learn  the  mechanics  of 
model  application.  The  first  example  is  wave  transformation  and  wave-current 
interaction  at  a  complex  coastal  inlet.  The  second  example  is  wind-wave 
generation  over  a  simple  bathymetiy.  The  first  example  focuses  on  wave 
propagation  and  the  second  example  on  wave  generation  (using  the  source  term 
options). 


Example  1 :  Wave  Propagation 
at  Ponce  de  Leon  Inlet 

The  first  example  is  based  on  a  field  site,  Ponce  de  Leon  Inlet,  located  on  the 
east  coast  of  Florida.  The  nearshore  bathymetry  at  Ponce  Inlet  is  complex 
because  of  the  presence  of  the  inlet,  a  navigation  channel,  jetties,  and  a  large  ebb 
shoal.  Tidal  currents  at  the  inlet  can  exceed  1  m/sec.  Engineering  problems  at 
Ponce  Inlet  that  require  information  about  the  nearshore  wave  field  include  scour 
in  the  inlet  throat  near  the  north  jetty,  erosion  of  the  north  spit  in  the  interior  of 
the  inlet,  migration  of  the  navigation  channel,  and  possible  breaching  shoreward 
of  the  north  jetty.  More  information  about  Ponce  Inlet  is  given  in  a 
reconnaissance  report  by  the  U.S.  Army  Corps  of  Engineers  (1993).  Also,  a 
field-data  collection  program  conducted  at  Ponce  Inlet  from  September  1995 
through  September  1997  is  documented  by  King  et  al.  (1999).  STWAVE 
applications  to  Ponce  Inlet  are  documented  in  Smith,  Militello,  and  Smith  (1998) 
and  Smith  and  Smith  (in  preparation). 


Model  input 

Model  parameters.  Figure  8  shows  the  bathymetiy  contours  within  the 
STWAVE  modeling  domain  selected  for  Ponce  Inlet.  The  offshore  boundary 
location  is  selected  at  an  approximate  depth  of  18  m.  This  depth  is  chosen 
because  the  bottom  contours  are  fairly  straight  and  parallel  at  this  depth  (out  of 
the  influence  of  the  ebb  shoal)  and  because  a  wave  gauge  (designated  as  DWGl 
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Figure  8.  Ponce  Inlet  bathymetry  and 
gauge  positions 


in  Figure  8)  was  located  near  this 
contour  to  provide  the  input  wave 
conditions.  The  lateral  boundaries  of 
the  domain  are  positioned  away  from 
the  influence  of  the  inlet,  to  areas  of 
fairly  straight  and  parallel  depth 
contours.  A  50-m  grid  cell  spacing  is 
selected  to  resolve  the  ebb  shoal  and 
inlet  bathymetry.  To  cover  the 
domain  with  50-m  resolution  requires 
53  cells  across  the  shore  (NI)  and 
1 12  cells  alongshore  (NJ).  Typical 
wave  periods  for  this  site  are  1 5  to 
5  sec  (0.0667  to  0.2  Hz).  To  resolve 
this  range,  30  frequency  bins  are  used 
with  an  initial  frequency  of  0.03 1  Hz 
and  a  frequency  increment  of 
0.0078  Hz  (range  of  frequencies  is 
0.03 1  to  0.258  Hz).  Because  the  wave 
propagation  distances  on  the  Ponce 
Inlet  grid  are  short  (2  km),  source 
terms  are  neglected  (IPRP  =  1 ).  Even 
with  strong  winds,  wave  growth  over 
the  short  distance  would  be  small. 
Wave-current  interaction,  though, 
could  be  significant  near  the  inlet,  so 
ICUR  =  1  to  include  wave-current 
interaction.  Selected  model  output 
locations  are  chosen  to  match  three 
directional  wave  gauge  positions  at 
Ponce  Inlet  (see  Figure  8).  The  depths 
at  the  three  gauges  are  approximately 
14,  7,  and  5  m.  The  gauges  are 


located  near  the  offshore  boundary,  on  the  ebb  shoal,  and  in  the  outer  inlet  throat 
(north  of  the  ebb  shoal).  The  parameter  input  file  for  this  case  is  given  in 
Appendix  A. 


Bathymetry.  The  bathymetry  contours  for  Ponce  Inlet  are  shown  in 
Figure  8.  The  grid  was  developed  using  a  grid  generator  (see  Chapter  5),  which 
accepts  random  x,  y,  and  depth  triplets  and  interpolates  them  onto  a  Cartesian 
grid  with  a  given  origin,  orientation,  and  resolution.  In  this  case,  high-resolution 
bathymetry  near  the  inlet  was  obtained  from  a  Scanning  Hydrographic 
Operational  Airborne  Lidar  System  (SHOALS)  (Lillycrop,  Parson,  and  Irish 
1996)  survey  and  supplemented  with  National  Oceanic  and  Atmospheric 
Administration  digital  bathymetry  (used  to  define  offshore  and  bay  bathymetry) 
and  shorelines.  The  north  jetty  at  Ponce  Inlet  is  represented  in  the  model  as  a 
series  of  land  cells.  The  Jetty  is  specified  as  being  three  grid  cells  wide,  which  is 
wider  than  the  real  structure.  Three  cells  are  required  to  block  wave  energy  from 
propagating  through  the  structure.  Future  upgrades  to  the  model  will  provide  the 
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capability  to  represent  structure  widths  that  are  less  than  the  grid  cell  spacing. 
The  bathymetric  input  file  for  Ponce  Inlet  is  given  in  Appendix  B. 

Incident  wave  spectrum.  The  incident  wave  spectrum  for  Ponce  Inlet  was 
generated  using  a  TMA  spectral  shape  (with  a  spectral  peakedness  parameter, 

Y  =  5.0),  cos""6  directional  distribution  (with  nn  =  12),  H„„  =  5.2  m,  =  12.8 
sec,  and  Q„,  =  -7  deg  (see  Table  1  for  guidance  on  spectra  shape  and  directional 
spreading  parameters).  The  input  spectrum  file  is  given  in  Appendix  C.  These 
incident  conditions  represent  extreme  high  wave  conditions  measured  at  Ponce 
Inlet  in  March  1 996. 

Current  field.  The  current  field  for  Ponce  Inlet  was  generated  using  the  tidal 
circulation  model  ADCIRC  (Luettich,  Westerink,  and  Scheffner  1992).  The 
current  field  output  from  ADCIRC  was  interpolated  onto  the  STWAVE  grid 
using  SMS  (see  Chapter  5).  The  input  current  field  is  given  in  Appendix  D.  The 
current  field  is  plotted  in  Figure  9. 


Results 

The  wave  parameters  for  the  three  gauge  positions  at  Ponce  Inlet  are  given  in 
the  selhts.out  file: 


Date 

1 

J 

^mo 

T-p 

e. 

96031202 

15 

21 

5.25 

12.8 

-6. 

96031202 

15 

65 

4.13 

12.8 

-5. 

96031202 

28 

49 

2.86 

12.8 

18. 

At  the  first  gauge  position  (I  =  15,  J  =  21),  the  wave  height  has  shoaled  slightly 
(from  5.2  m  at  the  boundary  to  5.25  m)  and  turned  slightly  more  shore  normal 
(from  -7  deg  at  the  boundary  to  -6  deg).  At  the  second  gauge  (I  =  1 5,  J  =  65), 
located  on  the  ebb  shoal,  the  wave  energy  has  dissipated  because  of  breaking 
(reduction  of  21  percent  in  wave  height  from  the  offshore  boundaiy).  At  the 
most  shoreward  gauge  (I  =  28,  J  =  49),  the  energy  has  been  dissipated 
significantly  (45-percent  reduction  in  wave  height  from  the  offshore  boundary) 
and  the  mean  direction  has  refracted  from  -7  deg  at  the  boundary  to  1 8  deg  in  the 
outer  throat.  The  wave  directions  at  this  gauge  are  turning  to  align  normal  to  the 
ebb  shoal  that  is  located  to  the  south  of  the  gauge  (positive  angles  are  more 
southerly  and  negative  angles  are  more  northerly). 

The  trends  of  wave-height  reduction  and  turning  of  the  wave  angle  can  also 
be  illustrated  by  examining  the  one-dimensional  wave  spectra  plotted  in 
Figure  1 0.  The  shape  of  the  spectra  stays  quite  similar,  but  the  energy  is  reduced 
because  of  depth-limited  breaking  from  the  offshore  to  the  outer  throat  gauge. 
Another  way  to  view  the  spectrum  is  to  integrate  over  all  frequencies  to  examine 
the  directional  distribution  of  the  wave  energy  (Figure  1 1 ).  The  directional 
distribution  narrows  as  well  as  reduces  in  energy  between  the  offshore  and  ebb 
shoal  gauges  because  of  refraction.  Between  the  offshore/ebb  shoal  gauges  and 
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Energy  Density  (m*m/rad)  <5  Energy  Density  (m*m/Hz) 


jre  1 0.  One-dimensional  wave  spectra  for  Ponce  Inlet 


Figure  1 1 .  Directional  distribution  of  Ponce  Inlet  wave  spectra 
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the  outer  throat  gauge,  the  mean  direction  shifts  by  25  deg,  again  because  of 
refraction. 

Figure  12  shows  a  contour  plot  of  the  wave  heights  over  the  entire  STWAVE 
domain.  The  gray  shade  contours  represent  wave  height,  and  the  background 
line  contours  represent  the  bathymetiy.  For  these  incident  wave  conditions,  the 
wave-height  variations  closely  follow  the  depth  contours,  because  depth-limited 
breaking  in  the  dominant  process  and  thus  the  wave  height  contours  mimic  the 
depth  contours.  For  cases  with  less  extreme  incident  wave  heights,  the 
maximum  wave  heights  are  found  on  the  ebb  shoal,  where  energy  is  focused  by 
refraction. 


Figure  12.  Wave-height  contours  for  Ponce  Inlet  example 
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Example  2:  Wind-Wave  Generation 

For  the  second  example,  consider  a  bay  that  is  rectangular  in  shape,  20  km 
east-west  and  10  km  north-south.  For  practical  purposes,  the  bay  is  isolated  from 
ocean  waves.  The  water  depth  is  a  constant  15  m,  with  a  gentle  side  slope  of 
1 : 100  at  the  eastern  end  of  the  lake  and  vertical  walls  on  the  other  three  sides. 
The  bay  has  a  rectangle  shape  with  simplified  bathymetry,  so  users  can  easily 
develop  the  input  bathymetry.  The  grid  domain  is  shown  in  Figure  13.  The  area 
of  interest  is  the  eastern  shoreline  of  the  bay,  where  a  project  is  planned  to  mine 
sand  to  construct  a  beach.  For  a  project  such  as  this,  wave  information  may  be 
required  to  evaluate  sediment-transport  rates  at  the  borrow  site  to  estimate 
infilling  rates  or  evaluate  longshore  and  cross-shore  sediment-transport  rates  at 
the  newly  constructed  beach  to  estimate  the  project’s  life  and  renourishment 
requirements. 


N 


Grid  Domain 


Borrow  Site 
(1=85,  J=25) 


Figure  1 3.  Schematic  of  Example  2  grid  domain 


Model  input 

Model  parameters.  Both  wave  propagation  and  source  terms  (IPRP  =  0)  are 
required  for  this  example,  but  currents  are  negligible  (ICUR  =  0).  The  special 
output  points  are  selected  to  examine  wave  growth  along  the  major  axis  of  the 
lake.  The  borrow  site  is  located  near  grid  cell  (85,25).  The  model  parameter 
input  file  is: 


0  0  10 
5  25 
25  25 
50  25 
75  25 
85  25 
95  25 

97  25 

98  25 
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99  25 

100  25 


Bathymetry.  The  bathymetry  and  shoreline  are  very  simple,  and  fine 
bathymetric  resolution  is  not  required.  A  grid  spacing  (DXINC)  of 200  m  is 
selected,  with  101  cells  in  the  cross-shore  (NI)  and  51  cells  in  the  alongshore 
(NJ).  To  simplify  the  input  file  for  this  example,  the  bay  is  assumed  to  have 
vertical  walls  on  its  the  north,  south,  and  west  shores  and  a  1:100  slope  on  the 
eastern  shore.  The  depth  in  all  grid  cells  along  the  entire  grid  boundary  is  set  to 
-1  for  land.  The  first  two  grid  rows  (J  =  5 1  and  J  =  50)  of  the  depth  file  are  given 
below: 


100  51  200. 
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15. 

15. 

15. 

15. 

15. 

15. 

13. 

11. 

9. 

7. 

5. 

3. 

1. 

-1. 

The  full  depth  file  would  repeat  the  second  series  above  for  rows  J  =  49  through 
2,  and  the  depths  for  row  J  =  1  would  again  all  be  -1.,  identical  to  the  first  series 
in  the  file.  Note  that  the  first  line  in  the  depth  file  represents  the  north  edge  of 
the  grid  (J  =  NJ),  reading  from  west  to  east  (I  =  1,NI),  and  the  last  series  in  the 
files  would  be  the  southern  edge  of  the  grid  (J  =  1). 

Incident  wave  spectrum.  Twenty  frequency  bins  are  used  with  an  initial 
frequency  of  0.04  Hz  and  a  frequency  increment  of  0.04  Hz.  This  coarse 
frequency  resolution  is  sufficient  because,  even  at  high  wind  speeds,  the  20-km 
fetch  will  limit  peak  wave  periods  to  approximately  5  to  6  sec.  The  input 
spectrum  is  veiy  simple  for  this  example  case.  Because  the  bay  is  treated  as  an 
enclosed  basin,  the  spectrum  on  the  “offshore”  boundary  (the  western  land 
boundary  for  this  case)  is  zero.  The  input  spectrum  includes  the  header  with  the 
date  and  wind  information,  followed  by  700  zeros  (20  frequencies  by  35 
directions).  For  a  date  identifier  of 980923,  a  wind  speed  of  20  m/sec,  a  wind 
direction  of  0  deg  relative  to  the  STWAVE  grid  (wind  blowing  from  the  west, 
for  this  example),  a  peak  frequency  of  0.8  Hz  (set  to  the  highest  frequency 
because  there  is  no  energy  in  the  input  spectrum),  and  0.0  m  water-elevation 
correction;  the  beginning  of  the  spectrum  file  is  given  by: 
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20  35 

0.04  0.08  0.12  0.16  0.20  0.24  0.28  0.32  0.36  0.40 

0.44  0.48  0.52  0.56  0.60  0.64  0.68  0.72  0.76  0.80 

980923  20.  0.  0.8  0. 

In  the  input  file,  these  lines  would  be  followed  by  700  zeros  (35  directions  by 
20  frequencies)  to  represent  an  incident  spectrum  with  no  energy. 

Current  field.  No  current  field  file  is  specified  for  this  case. 


Results 


The  results  of  greatest  interest  for  this  case  are  the  wave  parameters  at  the 
proposed  borrow  site  (I  =  85,  J  =  25).  The  wave  height  is  1 .52  m,  the  peak 
period  is  4.6  sec,  and  the  direction  is  0  deg  (from  the  west,  for  this  grid 
orientation).  The  output  file  with  selected  wave  parameters  is  given  below 
(header  has  been  added). 


Date 

1 

J 

^mo 

T-p 

0. 

980923 

25 

25 

0.86 

3.2 

0. 

980923 

50 

25 

1.22 

3.9 

0. 

980923 

75 

25 

1.45 

4.4 

0. 

980923 

85 

25 

1.52 

4.6 

0. 

980923 

95 

25 

1.55 

4.7 

0. 

980923 

96 

25 

1.53 

4.7 

0. 

980923 

97 

25 

1.50 

4.7 

0. 

980923 

98 

25 

1.47 

4.8 

0. 

980923 

99 

25 

1.23 

4.8 

0. 

980923 

100 

25 

0.59 

4.8 

0. 

These  output  wave  parameters  across  the  long  axis  of  the  bay  demonstrate  the 
wave  growth  as  a  function  of  distance  along  the  fetch.  These  wave  heights  and 
periods  are  plotted  as  a  function  of  distance  in  Figure  14.  The  wave  height 
grows  with  fetch  until  the  depth  decreases  on  the  eastern  shore  of  the  lake,  and 
the  wave  height  decays  because  of  depth-limited  breaking.  The  wave  direction 
down  the  long  axis  of  the  lake  is  aligned  with  the  wind  direction  (0  deg). 

Figure  1 5  shows  selected  one-dimensional  spectra  along  the  long  axes  of  the  lake 
(legend  provides  I,  J  locations).  Moving  from  west  to  east  (I  =  25  to  95),  the 
peak  frequencies  decrease  (peak  periods  increase)  and  the  total  energy  increases, 
as  the  waves  grow  from  input  wind  energy.  Between  I  =  95  and  1 00,  the  energy 
density  decreases  because  of  depth-limited  breaking.  These  one-dimensional 
spectra  were  calculated  from  the  spec .  out  file  by  integrating  the  spectra  over 
all  directions. 

To  evaluate  the  impacts  of  the  borrow  site  on  the  shoreline  or  the  sediment 
infilling  rate  at  the  borrow  site,  additional  STWAVE  runs  would  be  required. 
Typically,  the  local  wave  climate  would  be  examined,  and  representative  as  well 
as  maximum  wind  conditions  would  be  used  as  input  for  multiple  model  runs. 
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Figure  14.  Wave  height  and  period  as  a  function  of  distance  for  20  m/sec  wind 
speed 


Results  from  these  runs  would  then  be  used  statistically  to  evaluate  sediment- 
transport  potential  at  the  borrow  site  and  along  the  shoreline.  Note  that  the 
strongest  wind  speed  does  not  necessarily  lead  to  the  largest  waves  or  largest 
sediment-transport  rates.  Wind  direction  is  also  a  critical  parameter. 

This  example  of  wave  growth  was  for  a  completely  enclosed  basin.  The 
wave  growth  capabilities  of  the  model  are  also  applicable  on  the  open  coast.  For 
open-coast  applications,  the  input  spectrum  would  likely  not  equal  zero,  and  the 
offshore  and  lateral  grid  boundaries  would  not  be  land.  For  open  coastal 
applications,  the  local  wave  generation  (represented  by  the  source  terms)  is  often 
an  essential  process,  increasing  the  wave  height  near  shore  and  altering  wave 
directions. 
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Energy  Density  (m*m/H2) 


Figure  1 5.  One-dimensional  wave  spectra  for  20  m/sec  wind  speed 


■2 
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Appendix  A 

Sample  Model  Parameter  File 


The  parameters  required  to  run  STWAVE  are  specified  in  the  options .  std 
file.  A  description  of  the  file  and  formatting  is  given  in  Chapter  3  of  the  main 
text.  An  example  parameter  file  for  Ponce  de  Leon  Inlet,  Florida  (see  Example  1 
in  Chapter  6  -  main  text)  is  given  below: 

113 
15  21 
15  65 
28  49 


Appendix  A  Sample  Model  Parameter  File 


A1 


Appendix  B 

Sample  Bathymetry  File 


The  bathymetry  required  to  run  STWAVE  is  specified  in  the  depth  •  in  file. 
A  description  of  the  file  and  formatting  is  given  in  Chapter  3  of  the  main  text. 

An  example  of  the  first  seven  rows  (starting  with  J  =  NJ)  of  the  bathymetry  file 
for  Ponce  de  Leon  Inlet,  Florida  (see  Example  1  in  Chapter  6  -  main  text)  is 
given  below  (depths  in  meters  relative  to  mean  sea  level): 


53  112 
17.9093 
15.0553 
12.0976 
9.3289 
6.4268 
0.6400 
17.8912 
15.1025 
12.0504 
9.2965 
6.4597 
0.6400 
17.8730 
15.0376 
12-0228 
9.2506 
•6.2009 
0.6400 
17.8735 
14.9726 
11.8244 
9-3304 
6.2534 
0.6400 
17.9143 
14,9013 
11-7521 
9.1033 
6.0619 
0-6400 
17.9663 
14.8025 
11.6596 
8.9976 
6.0211 
0.6400 
17.9640 
14.8710 
11.5210 
9.0647 
5.7495 
0.6400 


50.0 

17.6657 

14.6698 

11.6441 

9.1748 

5.9270 

0.6400 

17.6474 

14.7452 

11.3595 

9.1355 

5.6139 

0.6400 

17.6292 

14.7832 

11.5448 

9.1210 

5.6831 

0.6400 

17.6320 

14.7182 

11.5183 

8.9854 

5.4282 

0.6400 

17.6802 

14.6412 

11.3491 

8.8021 

5.5081 

0.6400 

17.6959 

14.5424 

11.3985 

8.8395 

5.5151 

0.6400 

17.6947 

14.4960 

10.9378 

8.6888 

5.1866 

0.6400 


17.4214 

14.3101 

11.1662 

9.0352 

5.3055 

0.6400 

17.4032 

14.3855 

11.2053 

8.7766 

5.1999 

0-6400 

17.3849 

14.4608 

11.0923 

8.7173 

5.0017 

0.6400 

17.4084 

14.4639 

11.0874 

8.7941 

5.0103 

0.6400 

17.4390 

14.3811 

11.0092 

8.8039 

4.8664 

0.6400 

17.4407 

14.2823 

10.8634 

8.6251 

4.6782 

0.6400 

17.4410 

14-1835 

10.5820 

8.5083 

4.5709 

0.6400 


17.1925 

13.9923 

10.5809 

8.7189 

4.6770 

17.1589 

14.0510 

11.0999 

8.5754 

4.6779 

17.1414 

14.1097 

10.7910 

8.6526 

4.4674 

17.1786 

14.1765 

10.9693 

8-5154 

4.4061 

17.1826 

14-1086 

10.8806 

8.3929 

4.1354 

17.1867 

13-9758 

10.7462 

8.3809 

4.1240 

17.1872 

13.8431 

10.5504 

8.2821 

3.9764 


16.9141 

13.7952 

10.9137 

8.5227 

4.1037 

16.9146 

13.8539 

10.7788 

8.3478 

3.9818 

16.9177 
13. '9126 
10.6764 
8.4961 
3.7713 

16.9223 

13.8553 

10.3147 

8.3952 

3.5952 

16.9263 

13-7055 

10.6651 

8-2766 

3.4224 

16.9304 

13.5728 

10.5456 

8-0959 

3.2644 

16.9334 

13.4401 

10.3108 

8.1340 

2.9968 


16.6358 

13.5981 

10.5209 

8.1932 

3,3492 

16.6455 

13.6568 

10.4919 

8.2013 

3.2214 

16.6265 

13.6184 

10.4101 

8.1170 

2.9574 

16.6361 

13.4639 

10.3180 

8.0274 

2.8643 

16.6457 

13-3094 

10.4443 

7.9787 

2-6534 

16.6553 

13.1698 

10.2204 

8.0313 

2.5240 

16.6649 

13.0731 

10.1297 

7.8393 

2.8497 


16.3377 

13.3654 

10.4223 

7.9538 

2.7907 

16.3367 

13.3814 

10.4072 

7.9531 

2.6701 

16.2662 

13.2269 

10.3585 

7.7748 

2.4476 

16.2390 

13.0724 

10-2337 

7.7858 

2.2063 

16-2486 

12.9180 

10.0064 

7.5976 

2-1461 

16.2582 

12.7667 

9.9978 

7.6854 

2.0410 

16.2962 

12.7211 

10.0025 

7.5550 

2.0111 


16.0289 

13.0974 

10.2180 

7.7774 

2.2752 

16.0278 

13.0899 

10.1708 

7.5586 

2.1188 

15.9447 

12.9166 

9.9016 

7.5424 

1.9691 

15.8532 

12.7422 

10.0296 

7.3618 

1.7610 

15.8515 

12.5679 

9.8200 

7.3262 

1.6519 

15-8635 

12.3935 

9-9167 

7.3531 

1.5429 

15.9504 

12.2679 

9.6549 

7.1746 

1.4338 


15.7200 

15.4130 

12.7856 

12.4416 

10.0181 

9.5797 

7.3291 

6.9325 

1.6784 

1.0786 

15.7146 

15.3931 

12.7426 

12.3953 

9.7571 

9.4712 

7.2776 

6.7930 

1.5675 

1.0081 

15.6231 

15.3016 

12.6496 

12.3701 

9.8148 

9.4819 

6.9273 

6.8080 

1.3759 

0.8818 

15.5317 

15.2270 

12.4753 

12-2083 

9.6704 

9.5412 

7.0744 

6-5947 

1.2669 

0.7728 

15.4544 

15.1614 

12.3009 

11-8453 

9.5652 

9.4603 

6.9120 

6-6146 

1.1578 

0.6637 

15.5177 

15.1553 

11-8662 

11.2459 

9.4172 

9.3615 

6-8687 

6.6196 

1.0488 

0.6400 

15.6046 

15.2460 

12.0667 

11.4705 

9.3823 

9.2771 

6.7947 

6.3021 

0.9397 

0-6400 

Appendix  B  Sample  Bathymetry  File 


Appendix  C 
Sample  Incident  Wave 
Spectrum  File 


The  input  two-dimensional  spectrum  required  to  run  STWAVE  is  specified  in 
the  spec .  in  file.  A  description  of  the  file  and  formatting  is  given  in  Chapter  3 
of  the  main  text.  An  example  of  a  file  header  and  the  energy  densities  (meters 
squared/hertz/radian)  in  the  first  five  frequencies  (0.03 1 2  to  0.0624  Hz,  with 
35  direction  bins  per  frequency)  of  the  spectrum  for  Ponce  de  Leon  Inlet,  Florida 
(see  Example  1  in  Chapter  6  -  main  text)  are  given  below: 

30  35 

0.0312  0.0390  0.0468  0.0546  0.0624  0.0702  0.0780  0.0858  0.0936  0.1014 

0.1092  0.1170  0.1248  0.1326  0.1404  0.1482  0.1560  0.1638  0.1716  0.1794 

0.1872  0.1950  0.2028  0.2106  0.2184  0.2262  0.2340  0.2418  0.2496  0.2574 

96031202  0.  0.  7.81200E-02  0.340000 

0.00000  0.00000  0.00000  0.00000  0.00000  0.00000  0.00000 

0.00000  0.00000  0.00000  0.00000  0.00000  0.00000  0.00000 

0.00000  0.00000  0.00000  0.00000  0.00000  0.00000  0.00000 

0.00000  0.00000  0.00000  0.00000  0.00000  0.00000  0.00000 

0.00000  0.00000  0.00000  0.00000  0.00000  0.00000  0.00000 

0.00000  0.00000  0.00000  0.00000  0.00000  0.00000  0.00000 

0.00000  0.00000  0.00000  0.00000  0.00000  0.00000  0.00000 

0.00000  0.00000  0.00000  0.00000  0.00000  0.00000  0.00000 

0.00000  0.00000  0.00000  0.00000  0.00000  0.00000  0.00000 

0.00000  0.00000  0.00000  0.00000  0.00000  0.00000  0.00000 

0.00000  0.00000  0.00000  0.00000  0.00001  0.00003  0.00011 

0.00033  0.00082  0.00174  0.00326  0.00543  0.00810  0.01092 

0.01336  0.01488  0.01513  0.01404  0.01187  0.00912  0.00634 

0.00397  0.00222  0.00109  0.00046  0.00017  0.00005  0.00001 

0.00000  0.00000  0.00000  0.00000  0.00000  0.00000  0.00000 

0.00000  0.00000  0.00001  0.00006  0.00039  0.00182  0.00661 

0.01951  0.04837  0.10330  0.19343  0.32176  0.48009  0.64709 

0.79177  0.88232  0.89697  0.83206  0.70350  0.54069  0.37614 

0.23537  0.13135  0.06460  0.02756  0.00998  0.00297  0.00070 

0.00012  0.00001  0.00000  0.00000  0.00000  0.00000  0.00000 

0.00000  0.00000  0.00003  0.00036  0.00238  0.01123  0.04081 

0.12041  0.29856  0.63764  1.19397  1.98611  2.96345  3.99427 

4.88737  5.44632  5.53671  5.13607  4.34248  3.33753  2.32181 

1.45289  0.81075  0.39874  0.17010  0.06158  0.01835  0.00431 

0.00075  0.00009  0.00001  0.00000  0.00000  0.00000  0.00000 
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Appendix  D 

Sample  Current  Field  File 


The  input  current  field  required  to  run  STWAVE  is  specified  in  the 
current .  in  file.  A  description  of  the  file  and  formatting  is  given  in  Chap¬ 
ter  3  of  the  main  text.  An  example  of  the  first  row  (J  =  NJ)  of  the  currrent  file 
for  Ponce  de  Leon  Inlet,  Florida  (see  Example  1  in  Chapter  6  -  main  text)  is 
given  below  (u  and  v  velocity  components  in  meters/second): 


53  112  50.0 

96031202 
-5.9655029E-03 
6.6299401E-‘02 
-4-2588841E-03 
.  6.6767417E-02 
-1.7125718E-03 
6.4384326E-02 
1.6661584E-03 
6.2195450E-02 
3.2027178E-03 
6.0384035E-02 
4.2156633E-03 
5.8182545E-02 
7.3386896E-03 
5.4892469E-02 
8.3248559E-03 
4.6927348E-02 
7.1142977E-03 
3.5118129E-02 
5.0333366E-03 
2.5806203E-02 
O.OOOOOOOE+00 
O.OOOOOOOE+00 


6.6245958E-02 

-4.9415343E-03 

6.6381201E-02 

~3.1688120E-03 

6.5401770E-02 

4.7177821E-04 

6.2946498E-02 

2.5883764E-03 

6.1072327E-02 

3.8230885E-03 

5.9363544E-02 

6.0829669E-03 

5.6374207E-02 

8.4857941E-03 

5.0295375E-02 

7.6967701E-03 

4.0876567E-02 

5.8720950E-03 

2.8710742E-02 

4.2738337E-03 

O.OOOOOOOE+00 


-5.6242347E-03 

6.6327110E-02 

~4.0070452E-03 

6.6419885E-02 

-9.8446012E-04 

6.3874617E-02 

1.9737110E-03 

6.1821289E-02 

3.4305919E-03 

6.0042873E-02 

4.8182383E-03 

5.7583157E-02 

7.8893714E-03 

5.3458244E-02 

8.1154760E-03 

4.5243401E-02 

6.7089610E-03 

3.2525599E-02 

4.7800979E-03 

2.4354232E-02 

O-OOOOOOOE+00 


6.6273674E-02 

-4.6002306E-03 

6.6551067E-02 

-2.4407078E-03 

6.4892732E-02 

1.1998974E-03 

6.2570527E-02 

2.8955825E-03 

6.0723398E-02 

4.0193945E'-03 

5.8781888E-02 

6.7153089E-03 

5.5682164E-02 

8.5342079E-03 

4.8611302E-02 

7.4873865E-03 

3.7997276E-02 

5.3928681E-03 

2,7258173E-02 

4.0205559E-03 

O.OOOOOOOE+00 


-5.2828006E-03 

6.6354163E-02 

~3.7839543E-03 

6.5910846E-02 

-2.5634654E-04 

6.3366205E-02 

2.2812057E-03 

6.1447740E-02 

3.6268979E-03 

5.9703503E-02 

5.4505821E-03 

5.6983806E-02 

8.3813462E-03 

5.1979322E-02 

7.9061203E-03 

4.3559320E-02 

6.3036624E-03 

3.0454256E-02 

4.5269663E-03 

2.2901440E-02 

O.OOOOOOOE+00 
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